Local energy approach to the dynamic glass transition 
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We propose a new class of phenomenological models for dynamic glass transitions. The system 
consists of an ensemble of mesoscopic regions to which local energies are allocated. At each time 
step, a region is randomly chosen and a new local energy is drawn from a distribution that self- 
consistently depends on the global energy of the system. Then, the transition is accepted or not 
according to the Metropolis rule. Within this scheme, we model an energy threshold leading to a 
mode-coupling glass transition as in the p-spin model. The glassy dynamics is characterized by a 
two-step relaxation of the energy autocorrelation function. The aging scaling is fully determined 
by the evolution of the global energy and linear violations of the fluctuation dissipation relation are 
found for observables uncorrelated with the energies. Interestingly, our mean-field approach has a 
natural extension to finite dimension, that we briefly discuss. 
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Recent numerical studies have shed new light 

on the close connection between the dynamical slowing 
down of supercooled liquids and the topography of the 
underlying potential energy surface j^. In particular, a 
topological interpretation of the so-called mode-coupling 
temperature Tm has been confirmed H, Q. That is, 
for temperatures T smaller than Tm, the long-time relax- 
ation is governed by activated processes between basins 
of potential energy minima. For T > Tm , the thermal en- 
ergy is larger than the potential energy barriers, which al- 
lows the system to freely explore its configuration space. 

On the other hand, the (mean-field) mode-coupling ap- 
proximation for hquids 1 is exact for some disordered 
models p] such as the p-spin model 0- The spherical 
version of the latter, endowed with a Langevin dynam- 
ics, gives a clear illustration of the interplay between the 
dynamics and the topography of the energy landscape 
13 • The picture is the following 0|. The station- 
ary points whose energy is above the so-called threshold 
value Et are mainly saddles, i.e. the minimum value of 
the Hessian eigenvalues, Am, is negative. Those whose 
energy is below are mainly minima, i.e. Am > 0, and at 
Et the energy minima look flat: Am = 0. At low tem- 
perature, the expected equilibrium value of the energy is 
smaller than Et. Therefore, starting from high tempera- 
ture conditions, the system does not manage to reach the 
equilibrium because the phase space becomes flatter and 
flatter in the vicinity of Et- In the thermodynamic limit, 
the energy slowly drifts toward Ex but never reaches it, 
which leads to aging. The mode-coupling temperature 
corresponds to an equilibrium energy equal to Et- In 
a finite size system, the energy can cross the threshold 
so that the dynamics becomes activated. In this situa- 
tion, the phenomenological trap models 0, gi'^^ a 
fair description of the aging behavior. In particular, the 
regime just before the thermalization is expected to be 
well described jl^ by the activated trap model . 

The present Letter proposes a phenomenological mod- 
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FIG. 1: Schematic view of the local energies. A random Ising- 
like spin (f, J, ) is initially assigned to each site i. 



eling of the mode-coupling glass transition in the spirit 
of the spherical p-sm model. On one hand, we aim at 
providing new statistical tools to describe, at least at a 
mean-field level, the topology of an energy landscape, and 
more specifically the presence of an energy threshold. On 
the other hand, the present model has a natural extension 
to finite dimension, that we discuss in the conclusion. In 
this sense, we aim at giving a general framework to in- 
vestigate the passage from a mean-field description to a 
necessary real space description 

[3, El El of both the 

glass transition and glassy dynamics. 

Unlike the trap models that deal with global energies, 
our approach is based on the stochastic evolution of local 
energies {ei}i^i,,,N- These energies can be thought of 
as the energies of interaction felt by particles or spins, 
or as the energies of meso scop ic regions just as in the 
facilitated spin approaches [131 -see Fig. ^ For the sake 
of generality, we shall refer to the energy of a " site" . To 
further measure the response of neutral observables [IH j 
i.e. uncorrelated with the energies, a random Ising-like 
spin Si is initially assigned to each site i. Next, Si flips if 
and only if a transition is accepted for . 

The e'^s are supposed to evolve independently from 
each other. Thus, we drop the subindex i and deal with 
a single-energy e. In the spirit of the random energy 
model J^j, e{t) takes values in an infinite ensemble £ of 
independent random energies distributed according to a 
Gaussian distribution p(e) = exp(— e^) -where the energy 
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unit is set to one. Starting from an energy e at time t, the 
dynamics consists in drawing a new energy e' following a 
connectivity function /^(e — > e'), and then, accepting it 
or not according to the MetropoHs rule. The function 
depends on the mean energy per site e = ^le-i/N and 
reflects the topological properties of the energy surface 
at the global energy Ne. Roughly speaking, e can thus 
be thought of as an analog of the self-consistent magne- 
tization in the ferro-magnetism a la Curie- Weiss. 

We shall endow with a specific dependence on e to 
model the effect of an energy threshold. In this case, wc 
find at low temperature a spin/energy relaxation with 
two time-sectors (see Fig. OJ. The aging scaling is fully 
determined by the evolution of e and the spin observ- 
able gives linear fluctuation dissipation relations (FDR) 
that lead to effective temperatures larger than T. The 
energies have non-linear and non-monotonic FDR's. 

Mimicking a dynamic glass transition is con- 
strained by the following topological relation: 



p(e)/e(e ^ e') = p(e')/.(e' ^ e) 



(1) 



This is the reciprocity property of the connectivity. It 
says that the connection from a state A to a state B is 
identical to the connection from B to A. In the following, 
we consider a Gaussian/(5 function (see Fig. 01: 



((1 - e/edTny^e if e > e<j (2) 

6{e' + e) if e<ed (3) 



which verifies J^l if a = ■\/l — (1 — tjcdf . e^, an en- 
ergy that does not depend on the temperature and a 
positive exponent, are free parameters. 

To justify this choice, let us consider the Monte-Carlo 
evolution of some disordered spin model S where the 
spins are Ising-like. At each step, a spin is chosen with 
probability 1/iV and flipped or not according to the 
Metropolis rule. At infinite temperature, e = in our 
case, any transition is accepted. Then, between two suc- 
cessive flips of the same spin, the global energy, and hence 
the local energy, strongly fluctuate. We then expect a 
random choice for the local energies, which corresponds 
to /e=o(e e') ^ exp(— e'^) in our formalism. Close to a 
" mode-coupling" energy threshold, the scenario is differ- 
ent. For instance, in the p-spin model the closer to the 
threshold, the longer the system stays in the same region 
of the phase space 0. In other words, fast fluctuations 
do not drift the system away from an initial configura- 
tion. In the spin model S, fast fluctuations correspond 
to spin-flips and a configuration can be labeled by the 
data of its local energies. Therefore, close to the energy 
threshold, between two flips of the same spin, the cor- 
responding local energy, e, should not change much. A 
way to include this effect is to consider a trapping into 
a two-state process (TSP) with energies {e,— e}. This 
justifies that /^(e —>■ e') becomes more and more peaked 




FIG. 2: e-dependence of the connectivity function. For 
e ~ Bd, the evolution is essentially trapped within the TSP's 
{— e, e}. The greater the exponent v, the more narrow /e, i.e. 
the more efficient the trapping within the TSP's. e < in 
this figure. 



around — e as e — > Cd, Cd thus playing the role of the 
threshold -see Eq. Q and Fig. [3 At the energy thresh- 
old, the system consists of an ensemble of TSP's that are 
described by the (5- function ©. 

At high temperature, the Boltzmann equilibrium is 
recovered and e = — I/2T. Thus, after a temperature 
quench below Td — l/2|erf|, the energy first rapidly de- 
creases toward and then, is slowed down by the trap- 
ping of the TSP's. At low temperature, we numerically 
find that the time needed to reach (measured in num- 
ber of sweeps) is larger than N when i/ > 2, so that it 
diverges in the thermodynamic limit. 

Correlation functions We have numerically inves- 
tigated the connected energy/spin autocorrelation func- 
tions, Ce(t„,i) = E,((e»(i™)e^(t)) - (e,(t„))(e,(t)))/iV, 
and Csit^.t) = E^i{s^{t^)s^{t)) - {s^{t^)) (s,{t))) / N re- 
spectively. The brackets stand for an average over the 
noise history. The initial conditions are always taken at 
infinite temperature and the waiting time before mea- 
surements is noted t^. The simulations reported here 
were done with v = A although similar results are ob- 
tained for different v >2. 

Fig. O shows the spin autocorrelation at different tem- 
peratures. Three regimes must be distinguished. At high 
T (data not shown) , the system exponentially relaxes and 
the evolution is time translational invariant (TTI). At 
lower T >Td (Fig. |2K) , the autotocorrelation develops a 
plateau that becomes TTI at finite time (not diverging 
with N). At further lower T < Td, the system enters in 
the so-called aging regime (Fig. ^p). Then, we observe 
a short-time regime {t — tu, <ti tm) that becomes TTI 
as tw increases whereas for long times t — tu, 3> , the 
relaxation depends on tw 

These behaviors can be rationalized within a master 
equation (ME) approach. As e — > Cd, the distribution 
can be expanded in powers of (1 — tjed)^ . Keeping only 
the first terms, the ME for the energy e reads 20] : 

dtP{e,t) = ~w{e ^ -e)P{e,t) + w{-e ^ e)P{-e,t) 
-(I - e(t)/edrF, (P(±e), P'(-e), P"(-e))(4) 

with e{t) — J deeP{e,t). Fi{-) is a functional of P{x,t) 
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FIG. 3; Spin autocorrelation, Cd = -0-7 ^ Td ~ 0.71. (a) 
Equilibrium curves. N = 2 x 10^. Inset: the time is rescaled 
by (1 - T/Tdy-, the curves are normahzed at the plateau 
value, (b) aging behavior. T = 0.6. A'^ = 5000. Dashed 
lines are the plateau values predicted for a full decoupling of 
the timescales (see text). Similar curves are obtained for the 
energies. The times are measured in number of sweeps. 



and w is the Metropolis rate. 

When dealing with the above two-time autocorrelation 
functions, the initial conditions in Eq. must be taken 
at time i^. Thus, the two first terms of the right hand 
side, that account for the relaxation within the TSP's, 
contribute to the fast relaxation toward the plateau. The 
last term allows the system to further relax with a typical 
timescale ^ (1 — e{t)/ed)~'' that can be time-dependent 
if t{t) does not reach a stationary regime. Therefore, at 
equilibrium the relaxation time for the decay from the 
plateau should scale like (1 — e/ed)~^ ~ (T — Td)~^ as 
T ^ Td since the closer to Cd, the sharper the timescale 
separation in Q). Fig. |2t confirms this scaling law once 
normalized the plateau values. 

For T < Td, e{t) drifts toward without reaching it. 
In the limit — > oo (rigorously taken after N — > oo), fast 
and slow timescales totally decouple. First, the TSP's 
locally equilibrate, which leads to P{e,t)w{e — * — e) ~ 
P{—e,t)w{—e — » e). Next, under these conditions the 
ME (gl) reduces to: 

dtP{e, t) = -(1 - e{t)/edrF2{P{e),P'{e),P"{e)) (5) 

where F2 is a functional different than Fi |]20| . This gives 
the dynamical evolution of -P(e, t) during the decay from 
the plateau. The timescale (1 — e{t)/ed)^'^ is now time- 
dependent. Furthermore, considering the decoupling of 
fast and slow timescales, the spin autocorrelation func- 
tion in the aging regime can be written as [20j : 



Cs{tw,t) = j dewdeG{e,t\en,,tn,){p{e) - p{-e)) {Q) 

X {p{ew) - p{~e.^))G{e^,t^) 

In this relation, G(e, t) is the energy density of the TSP's 
and is given by {P{e,t) + P{—e,t))/2. G{e,t\ew,tw) is 




FIG. 4: Energy autocorrelation vs. 4>t^{t). N — 4000. 
Ed ~ —0.8 ^ Td = 0.625. For each simulation, t runs from tw 
to 1000 X tjn. The curves are normalized at the plateau value 
with the shortest t^. (1 — t{t)/ed)~'^ was fitted by stretched 
exponentials to obtain 4>t^{t). Exponents of the stretched 
exponentials are pretty similar and are around 0.1. The fit 
correlation coefficient was always above 0.9995. Similar re- 
sults are obtained for the spin autocorrelation. 



the corresponding propagator with the initial condition 
P(±e,t) = S{e =F enj)S{t — tw). p[e) is the Boltzmann 
weight restricted to the TSP {e, — e}. Thus, the relation 
© means that the decay from the plateau comes from 
the decorrelation of the TSP's {e, — e} due to a diffusion 
of the energy e governed by the equation |SJ. 

Denoting (a;) = /^^ dt' {l~e{t') / edy , Eq. © calls 

for a solution G{e,t\ew,tw) = G(e,(/)t^(t)|etu,0). Insert- 
ing this relation into Q and noticing that G{e,tw) be- 
comes stationary as ^ 00, the aging scaling reads: 



(7) 



the reasoning for the energies being identical. This scal- 
ing is numerically well verified (see Fig. which cor- 
roborates our treatment of the aging regime. Moreover, 
in this scope, the plateau value for the spins at time tw is 
equal to / de G{e, tw){p(e) —p{—e))^, that is also numer- 
ically well verified (see Fig. |SJ). The Ed wards- Anderson 
parameter qea is given by this quantity when tw ^ 00. 

To summarize, two distinct timescales appear at low 
temperature. A fast one coming from the relaxation 
within the TSP's, and a slower one that increases as e(t) 
drifts to Cd. Notice that a full aging regime {t/tw scaling) 
is expected only if (1 — e{t)/ed)~'' decreases as 1/t since 
then (f)t^ (X \og(t/tw). However, we have never seen this 
regime in our simulations. 

Response to a field The two-time magnetic suscep- 
tibility corresponding to the spin autocorrelation func- 
tion reads Xs{tw,t) = d{^s{t)) /dh{tw)\h=o- {^s{t)) = 
jfJ^i^iSiit)) is measured a time t — tw after having 
switched on, at time tw, a small constant field h cou- 
pled to ^s. The fi's are quenched random variables that 
take values ±1 |23|. At low temperature and in the long- 
time limit, a piecewise relation between Xs and Cs has 
been found in the p-spin model 



Xsitw,t) = il-Csitw,t))/T if Gs>qEA 



(8) 
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FIG. 5: FDR. iV = 
equilibrium relation, 
ing h = 0.1. 
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1000, Ed = -0.7. The solid line is the 
The responses have been computed us- 



- iqEA-C'sit^,t))/T, if Cs<qEA (9) 

The relation JHJ is the usual equilibrium relation and is 
due to a local equilibrium property of the system. The 
linear relation Q in the aging regime leads to an effec- 
tive temperature Tg >T that shares the thermodynamic 
properties of canonical temperatures 2M ■ Interest- 
ingly, such a behavior has been later observed in numer- 
ical simulations and in experiments as well '2l'|, which 
prevents any artifact of the mean- field treatment. 

We computed Xs in our system by adding to a mag- 
netic energy e,„ — —h^iSi. Fig. [3 shows, for different 
t^, the plot Tx(tu),t) — u{C{tyj,t)) obtained in the lin- 
ear regime of the response. We see that it can be di- 
vided into two sectors corresponding to the fast relax- 
ation and the aging regime respectively. The former gives 
an equilibrium-like relation. This directly results from 
the local equilibrium within the TSP's (see above). The 
interesting results come from the aging part since it seems 
that we have a linear FDR for which we can define an 
effective temperature Tf. > T. As in the p-spin models, 
Te slightly increases as T decreases. 

We have also investigated the response of the local en- 
ergies. We found non-monotonic behaviors in the aging 
regime similarly to the results of in finite di- 
mension kinetic constrained models. Thus, our results 
suggest that even in mean- field situations, local energies 
may respond non-monotonically which may lead in some 
case to negative Te's [25j |. 

Conclusion We have considered the stochastic dy- 
namics of N independent local energies. At each Monte- 
Carlo step, the new energies are drawn from a distribu- 
tion /<; whose properties self- consistently depend on the 
global energy Nt of the system. The dependence is cho- 
sen in order to model an energy threshold in the spirit of 
the p-spin model. Given this mode-coupling scenario, we 
obtain at low temperature the typical properties of ag- 
ing. In particular, our approach gives a rather intuitive 



(and simple) illustration of how two time-sectors may ap- 
pear in a mean-field situation. The linear violation of the 
fluctuation dissipation relation, for observables uncorre- 
lated with the energies, further confirms the similarity 
between our phenomenological description and the mi- 
croscopic disordered models. 



Interestingly, one can generalize this mean-field ap- 
proach to finite dimension replacing e by a local energy 

~ S<j> ^0 where the sum is taken over the z near- 
est neighbors. In this case, one expects Td to become a 
crossover below which dynamical heterogeneities should 
play an important role. 
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